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Asymmetric exclusion processes with locally reversible kinetic constraints are introduced to in- 
vestigate the effect of non-conservative driving forces in athermal systems. At high density they 
generally exhibit rheological-like behavior, negative differential resistance, two-step structural relax- 
ation, dynamical heterogeneity and, possibly, a jamming transition driven by the external field. 



There is a growing appreciation that glassy relaxation 
can be ascribed to purely dynamic restriction on the par- 
ticle motion with static correlations and related thermo- 
dynamic factors playing little or no role [H, @, Q ■ Kinet- 
ically constrained models (KCMs) provide a simple and 
elegant way to rationalize this idea and have spurred the 
challenge to reproduce much of the observed glassy be- 
havior [J], including detailed predictions that first origi- 
nated in apparently unrelated and complementary mean- 
field approaches [H, While investigations of KCMs 
have been mostly focused on equilibrium and aging dy- 
namics, there are relatively few studies dealing with the 
effects of nonconservative forces. The issue is of spe- 
cial interest in rheology, where the apparent viscosity of 
structured liquids is found to depend on the applied stir- 
ring force [7(. Shear-thinning refers to the commonly 
observed situations in which the viscosity decreases at 
increasing forces, and is well described by mode-coupling 
theory [8|. The opposite, shear-thickening, behavior is 
less common and much more difficult to predict. In 
some cases the viscosity increase can be so dramatic 
that the macroscopic motion may even stop, the liquid 
jams @, 0]- It has been suggested that shear-thickening 
and jamming (loj are related to an underlying glass tran- 
sition [ll|. The idea has been interpreted microscopi- 
cally in terms of an entropy-driven inverse freezing [12j | 
and recent experiments have come to support this in- 
terpretation 3, HI- The latter approach, however, re- 
quires some special thermodynamic and structural fea- 
tures which are absent in athermal shcar-thickening sys- 
tems (e.g. concentrated suspensions of hard spheres). 
Accordingly, any attempt at unifying different types of 
dynamical arrest in systems dominated by steric hin- 
drance and packing effects should explain the simultane- 
ous emergence of thinning and thickening behavior with 
their underlying glassy dynamics. 

In this Letter I will show that rheological-like behavior 
occurs in microscopic models of finite dimensional par- 
ticle systems interacting only through non-Hamiltonian 
forces and purely kinetic constraints. Evidence is pro- 
vided by considering a variant of the asymmetric simple 
exclusion process (ASEP) [l^|, in which particle motion 
obeys an additional constraint motivated by lattice glass 
models [l6| . Similarly to the ASEP the transition proba- 
bilities satisfy local detailed balance although, due to the 



periodic boundary condition, the driving field cannot be 
derived from a Hamiltonian. The applied field induces a 
nonequilibrium steady state (NESS) in the system. At 
small field the transport is ohmic and thinning behavior 
is observed. At high density, due to the imposed kinetic 
restrictions on the particle motion, the ohmic regime 
shrinks and the current decreases at increasing field, a 
negative differential resistance effect. The system can 
thus be driven, at constant density, into a regime close 
to jamming in which the dynamics is slow but station- 
ary. In this regime, shear-thickening and other salient 
features of glassy behavior like two-step relaxation and 
dynamical heterogeneities are observed. 

The model. ASEP describes a ID lattice model 
of hard-core particles hopping randomly to a vacant 
nearest-neighbour site with rates p± depending on the di- 
rection of the particle move. Current-carrying NESSs are 
maintained either by fixing periodic boundary conditions, 
or by connecting the open boundaries to two particle- 
reservoirs at different densities. In the latter form the 
ASEP was originally introduced to model the biopoly- 
merization kinetics on nucleic acid templates [17j . and 
has come to play a paradigmatic role in recent develop- 
ments of nonequilibrium statistical mechanics [lil 
ASEP-type models have been also studied in relation 
to vehicular traffic, molecular motors and intra/extra- 
ccllular transport [l9(. 

In this work, a variant of the 2D ASEP with discrete- 
time evolution is introduced in which particle hopping 
occurs with probability p — min{l, exp (E ■ dr)}, where 
E is the applied field and dr is the displacement unit 
vector, provided that an additional condition is met: the 
particle is constrained to have a number of neighboring 
particles below a certain threshold both before and af- 
ter the move (ltjj . The constraint is naturally inspired 
by the cage effect in viscous liquids and leads to glassy 
relaxation at high density. In zero field, one recovers the 
boundary-driven constrained diffusion model of Ref. [20| . 
For periodic boundary conditions, the kinetic nature of 
the constraint guarantees that the equilibrium measure 
is trivial. The fully irreversible case is obtained in the 
opposite limit of infinite drive, corresponding to the ki- 
netically constrained version of the totally asymmetric 
simple exclusion process. 

For simplicity we consider square-lattice systems (of 
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size L 2 ) with periodic boundary conditions in which the 
driving force is applied along a lattice axis. The kinetic 
constraint is chosen to depend isotropically on the nearest 
neighbors with a threshold set to three. This constraint 
ensures that diffusion coefficient at equilibrium is finite at 
any nonzero vacancy [![. Thus, no "true" glass transition 
is present in zero field although the dynamics is slow at 
high density and characterized by the cooperative motion 
of spatially extended mobility defects, whose size grows 
with the density [l[ . For finite field, every allowed parti- 
cle move in the direction opposite to the field can always 
occur with finite probability. This would imply that the 
irreducibility and ergodicity proofs given in Ref. [H may 
be extended to the NESS, provided that the field inten- 
sity E is finite. Some technical steps, however, may be 
not straightforward because, in contrast with the equi- 
librium case and with the ASEP, the probability weight 
of constant-density configurations is nonuniform in the 
NESS, even though the average density profile is flat. 

A first insight in the flow behavior can be gained by 
looking at the particle current J vs the driving force. 
At densities below p ~ 0.79, J first grows linearly with 
E (ohmic regime) and then tends to saturate at large 
fields. This behavior is qualitatively similar to that ob- 
served in the ASEP and will not be further considered 
here. At densities above p ~ 0.79 the transport proper- 
ties differ rather markedly from the ASEP: the current 
shows a transition from the ohmic regime at small fields 
to a non-monotonic regime in which the current attains 
a maximum at finite field (rather than at saturation as 
in the previous case), and then decreases for larger field, 
see fig. f . The region with negative slope in the "current- 
voltage" characteristics is traditionally known as incre- 
mental negative resistance (NR) and is a key ingredient in 
many biological systems and solid-state devices. Micro- 
scopic stochastic models yielding NR are known pol 21 1 . 



In our case, NR does not depend on any static interaction 
and occurs because, at high density and increasing field, 
the particle rearrangements needed to remove obstruc- 
tion to the flow, require more and more particle moves 
against, or normal to, the field direction (see, fig. 4c in 
Ref. [22j for an illustration). Three distinct behaviors 
may arise at high p: 

I. Saturation current attains a finite (nonzero) value. 
II. Saturation current vanishes. 
III. Current vanishes at a finite driving force. 

Numerical results suggest that regime I occurs in the 
range 0.79 < p < 0.83 while regime II appears at higher 
density. For the present model, the ergodicity result men- 
tioned above should prevent the existence of regime III in 
the thermodynamic limit, for large but finite E. Distin- 
guishing between the two latter regimes, however, may 
be difficult due to the strong finite-size effects related 




FIG. 1: Stationary current J vs field E at density p. 

to bootstrap percolation [l], 0] ■ The characterization of 
these effects is notoriously difficult and will not be at- 
tempted here. Kinetic constraints having a more involved 
dependence on the neighboring particles mm, wm ar- 
guably have a jamming transition at finite field in the 
thermodynamic limit. Moreover, KCMs having single 
mobility defects may only have type-I behavior, no mat- 
ter how large the particle density. This is confirmed by 
the behavior of a noncooperative driven KA model on 
triangular lattice. 

NESS dynamics is expected to be slow in the NR region 
of high density as the current becomes vanishingly small 
at increasing forcing. This can be analyzed through the 
persistence function (j)(t), i.e. the probability that the 
occupation variable of a lattice site has never changed 
between times and t. The two stages responsible for 
equilibrium relaxation - which consist of particle rear- 
rangements within spatially extended mobility defects 
and their subsequent coalescence [l[ - are found to re- 
spond differently to the applied field. 

At small field, relaxation under flow is slower than 
at equilibrium on short and intermediate time scales. 
At longer time scales, on the other hand, the macro- 
defects coalescence speeds up, and the relaxation behav- 
ior changes smoothly from a stretched to a simple expo- 
nential decay. The two stages of relaxation are well sep- 
arated by a density-dependent crossover time, see fig. [2] 
(upper panel). 

At large field, the thinning regime disappears because 
the time over which the first-stage relaxation occurs in 
the NESS exceeds the one required for the second-stage 
to take place in equilibrium. Consequently, <f> develops 
a more marked two-step behavior in which the corre- 
lation decay becomes slower and slower at increasing 
field. The resulting thickening-like behavior is shown 
in the lower panel of fig. [2j The late stage relax- 
ation in this regime obeys the superposition principle, 
<j){t) = $(i/r rol (p, E)), where the scaling function $ is 
a simple exponential for the model under consideration, 
and r rcl (p,E) = J <f>(t)dt is the integrated relaxation 
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p = 0.86 




FIG. 2: Steady state persistence <j> vs time t at density p, for 
small (upper panel) and large (lower panel) applied field E. 



time. T te i is a measure of the system viscosity and en- 
codes the overall relaxation resulting from the superposi- 
tion of the two above-mentioned competing effects. Fig. [3] 
clearly shows that at large enough density, r rol first de- 
creases at small field and then increases (exponentially) 
at large field. To better appreciate the result above one 
should consider that outside the linear response regime 
the "resistivity" E/J increases with E, while the relax- 
ation time is roughly independent of E at low density. 
This confirms that viscosity (rather than "resistivity") 
does generally account for the genuine flow properties of 
our model. The nonmonotonic dependence of viscosity 
at high density and the related transition from thinning 
to thickening behavior, are features observed in sheared 
hard-sphere suspensions Q. In fact, flow curve analogs, 
obtained by identifying the shear stress with E and the 
shear rate with E/r lel {p, E), arc similar to those found in 
the "Model II" of the schematic mode-coupling theories 
studied in [ill ]. In this case, since kinetic constraints are 
"hard", no fluidisation occurs after jamming no matter 
how strong the force may be. If violation of constraints 
is allowed at large field, however, one should recover the 
flow curves of "Model I" and "Model III" of Ref. [ll|. 

Notice that in contrast with the nonmonotonic behav- 
ior of viscosity and current, the fraction of blocked par- 
ticles increases with E, even when transport is ohmic. 
The driving force thus generally enhances the clustering 
of particles and appears to be akin to a static short-range 
attraction. In fact, models of this type do have a non- 
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FIG. 3: Integrated relaxation time r Te \ vs driving force E. 



monotonic dependence of the relaxation time with the 
attraction strength at equilibrium 



25[. This observa- 



tion suggests an analogy between attraction-driven in- 
verse freezing and shear-thickening/jamming. 

Finally, dynamic heterogeneity in the NESS can be 
quantified by the fluctuations of persistence, Xi(t) = 
N(((t> 2 (t)) - (<P(t)) 2 ) where N = pVjM,\¥$. Similarly 
to equilibrium supercooled liquids [27( we find that dy- 
namic heterogeneity and long range order grows with the 
applied field, as exemplified by the increasing peak in 
the dynamic susceptibility X4(*)i an d consistently with 
the monotonic increase of the fraction of blocked particle 
with E. However, in the thinning regime the peak occurs 
at shorter and shorter times as the field grows. In the 
thickening regime, on the other hand, the peak height 
and the time at which the peak occurs both increases 
with the applied field, as shown in fig. 01 Interestingly, 
even though the dynamics is obviously nonisotropic, the 
longitudinal and transverse persistence (and related fluc- 
tuations) are almost indistinguishable. Tiny differences 
only appear in the early stage of relaxation at large field. 
Thus, there is no substantial change in the corresponding 
longitudinal and transverse relaxation times. 

Conclusions To summarize, when KCMs are driven 
into a nonequilibrium stationary state they show both 
an initial speed-up of the dynamics and then a pro- 
nounced slowdown at increasing fields. This suggests a 
common kinetic mechanism for the nonmonotonic vis- 
cosity of athermal systems and their thinning-thickening 
transition. The basic ingredient is the presence of two re- 
laxation stages responding differently to the applied field. 
Such a behavior occurs in the simplest case of a spatially 
uniform nonconservative force. We expect that including 
a space dependent driving (to mimick more realistically a 
shear stress) yields similar results. In particular, the pe- 
culiar time evolution of dynamic heterogeneities observed 
here in the thinning and thickening regime should be ex- 
perimentally accessible in sheared granular materials and 
concentrated colloidal suspensions. Finally, the existence 
of a jamming/blocking transition, possibly induced by 
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FIG. 4: Steady state four-point correlation Xi vs time t for 
small (upper panel) and large (lower panel) applied field E. 



strong finite-size effects, implies that driven cooperative 
KCMs, are able to sustain indefinitely an applied stress 
above a certain threshold, and thus to provide a micro- 
scopic realization of fragile matter 28|. 
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